Міністерство освіти і науки України Національний університет “Львівська політехніка” Інститут Комп’ютерних наук та Інформаційних технологій Кафедра ПЗ Звіт з лабораторної роботи №3 На тему: “ Формування вхідних потоків і контролю їх якості для комп’ютерного моделювання складних систем статистичними методами” з дисципліни “ Моделювання ПЗ ” Мета роботи: одержання послідовностей псевдовипадкових квазірівномірно розподілених чисел запропонованими методами та перевірка їх якості. Завдання: Одержати послідовності псевдовипадкових квазірівномірно розподілених чисел при заданих параметрах для кожного з нижчевказаних способів: методу серединних квадратів; мультиплікативного методу; за допомогою вбудованих функцій із програмного забезпечення комп’ютера. Кожну з генерованих послідовностей перевірити на рівномірність двома методами: за гістограмою; за непрямими ознаками. Кожну з генерованих послідовностей провірити на незалежність методом, основаним на обчисленні кореляційного моменту. Перевірити одержані послідовності на стохастичність. Код програми Метод серединних квадратів using System; using System.Collections.Generic; using System.ComponentModel; using System.Data; using System.Drawing; using System.Linq; using System.Text; using System.Windows.Forms; namespace Lab1 { public partial class Form1 : Form { public Form1() { InitializeComponent(); } private static int t = 1; private static int count = 20; #region Метод серединних квадратів private void button1_Click(object sender, EventArgs e) { List<double> mas = new List<double>(); try { Output1.Clear(); double X1 = double.Parse(textBox1.Text); int N = int.Parse(textBox2.Text); int length = X1.ToString().Length - 2; if (length % 2 == 1) length++; StringBuilder sb = new StringBuilder(); string format = "{0:F" + length * 2 + "} "; string format_small = "{0:F" + length + "} "; sb.AppendFormat(format_small, X1); double Xi = X1; for (int i = 1; i < N; i++) { try { string tmp; Xi = Xi * Xi; tmp = string.Format(format, Xi); Xi = double.Parse(" 0," + tmp.Substring(length / 2 + 2, length)); mas.Add(Xi); sb.AppendFormat(format_small, Xi); } catch { } } Output1.Text = sb.ToString(); } catch (System.FormatException exc) { MessageBox.Show("Помилка при введенні даних для методу серединних квадратів"); } #region обчислення кореляційного моменту int NN = mas.Count; double SXi = 0, SXit = 0, SXiXit = 0, SXi2 = 0, SXit2 = 0; for (int i = 0; i < NN - t; i++) { SXi += mas[i]; SXit += mas[i + t]; SXiXit += mas[i] * mas[i + t]; SXi2 += mas[i] * mas[i]; SXit2 += mas[i + t] * mas[i + t]; } double DXi, DXit; DXi = 1.0 / (NN - t) * SXi2 - 1.0 / ((NN - t) * (NN - t)) * SXi * SXi; DXit = 1.0 / (NN - t) * SXit2 - 1.0 / ((NN - t) * (NN - t)) * SXit * SXit; double Pt; Pt = (1.0 / (NN - t) * SXiXit - (1.0) / ((NN - t) * (NN - t)) * SXi * SXit) / Math.Sqrt(DXi * DXit); textBox3.Text = Pt.ToString(); double B = 1.0 / NN; ...